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Using the Runge-Gross theorem that establishes the foundation of Time-dependent Density Func- 
tional Theory (TDDFT) we prove that for a given electronic Hamiltonian, choice of initial state, and 
choice of fragmentation, there is a unique single-particle potential (dubbed time- dependent partition 
potential) which, when added to each of the pre-selected fragment potentials, forces the fragment 
densities to evolve in such a way that their sum equals the exact molecular density at all times. 
This uniqueness theorem suggests new ways of computing time-dependent properties of electronic 
systems via fragment-TD'DFT calculations. We derive a formally exact relationship between the 
partition potential and the total density, and illustrate our approach on a simple model system for 
binary fragmentation in a laser field. 



Time-dependent density functional theory (TDDFT) 
[U [2] allows one to predict, in principle, the evolution 
of the non-relativistic density n(rt) of a system of in- 
teracting electrons subject to a time-dependent external 
potential v(rt). Given an initial wave function, the time- 
dependent electron density determines the external po- 
tential up to a time-dependent constant (Runge-Gross 
theorem [RG] pQ) and the density may be found by solv- 
ing the time-dependent Kohn-Sham (TDKS) equations. 
These equations make it possible to perform practical 
calculations to propagate the electronic density and its 
related quantities such as linear response functions. Due 
to its wide range of applications, TDDFT is expected 
to continue being a workhorse in the coming years for 
chemistry, physics, and materials engineering [3]. 

Although the computational cost of TDKS calcula- 
tions is low compared to that of other many-body tech- 
niques, new ideas are needed to enable the study of 
larger systems with improved efficiency and accuracy. 
For the ground-state problem, 'divide-and-conquer' frag- 
mentation techniques have been developed [3] and ap- 
plied successfully through the use of readily-available par- 
allel computers. Related strategies have also been de- 
veloped recently for the time-dependent problem within 
TDDFT [5]. For example, Casida and Wesolowski [6 
introduced a methodology to perform time-dependent 
calculations within frozen-density embedding theory. It 
has been shown that this method yields better results 
than "supermolecular" techniques in some cases [7J [3]. 
Other extensions include linear-response TDDFT for 
molecules in solvents [§] and TDDFT for interacting chro- 
mophores |10j . Also, time-dependent calculations within 
subsystem-DFT have been reported (see [11] and refer- 
ences therein). In subsystem-DFT, the density of the 
system is split into densities of localized subsystems. 
Then the Kohn-Sham kinetic energy of the total system 
is truncated, and the energy functional is approximated 
by a functional of the subsystem densities; its minimiza- 
tion leads to Kohn-Sham equations for each subsystem. 
Neugebauer formulated this theory within linear response 



in the frequency domain and showed that it yields results 
consistent with conventional TDDFT [T2]. For the case 
of dissipative dynamics, Zhou et al. |13) showed how the 
RG theorem can be applied and Kohn-Sham equations 
developed for open systems given an initial state, mem- 
ory kernel, and system-bath correlation. 

Among density-based ground-state fragmentation 
techniques, Partition Density Functional Theory (PDFT) 
[T3] is a reformulation of Density Functional Theory that 
allows one to find the solution to the KS equations with- 
out solving the total molecular problem directly. The 
idea is to partition the external potential into an arbi- 
trary number of fragment potentials. The total energy 
of the isolated systems is minimized under the constraint 
that the fragment densities sum to the correct molecu- 
lar density. The Lagrange multiplier associated with the 
constraint (i.e. the partition potential) can be found by 
inversion if the total density is known |15j . or via the 
self-consistent procedure of Ref. [T3] if it is not. Ev- 
ery fragment is subject to the same partition potential. 
In contrast with quantum mechanical embedding theo- 
ries (except for the latest version of quantum embedding 
[IB] ) and with subsystem-DFT, this potential is global 
and unique [17]. The set of fragment densities obtained 
for a given choice of external-potential partitioning is also 
unique. As Pavanello [TH] recently suggested, this unique- 
ness feature of PDFT makes it a suitable candidate to 
simplify the formulation of subsystem-DFT. This letter 
reports on foundational work for such developments. We 
extend PDFT to the time-dependent regime and show 
how the time-dependent external field can be partitioned. 
A new potential termed the time- dependent partition po- 
tential is introduced in the formalism in order to repre- 
sent the exact time-dependent electronic density. 

To extend PDFT to the time-dependent domain we re- 
call that there is no minimum principle from which the 
TDKS equations can be derived [HI [20]. In view of this, 
we follow a deductive approach to define our TDKS equa- 
tions. Our goal is to provide a fragment-based solution to 
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the Liouville equation (we use atomic units throughout) 



■ dT(t) 
' dt 



H v (t),T{t) 



(1) 



If r is a pure density matrix then Eq. (|1J) is equiva- 
lent to the time-dependent Schrodinger equation. We 
suppose that the initial state T(t ) is given. In stan- 
dard DFT notation, the Hamiltonian is given by H v (i) = 
T + Vee + J d 3 r v (rt) h (r). It is convenient to express the 
external potential v(rt) as the sum of the potential v(r) 
due to the M nuclei v(r) = — J2 a ^a/\ v — R-a|, which is 
not explicitly time-dependent, and an additional poten- 
tial VE^t) containing all of the explicit time-dependence 
due to external fields: 



v(rt) = v(r) + VE(rt) 



(2) 



Our task is to divide the quantum system into Nf frag- 
ments of interacting electrons. This is done by assigning 
an external potential v a (rt), Hamiltonian H a (t), and ini- 
tial state T a (to) to each fragment. Out of the infinitely 
many ways to choose the fragment potentials, there are 
at least two cases that are physically relevant: (i) Direct 
partitioning of the time-dependent external potential in 
analogy to ground-state PDFT: v(rt) — Yla* v a( r t)- For 
example, if Nf = M, there are cases of interest where we 
could define v a (rt) = —Z a /\r — Tl a (t)\. In such cases, 
the electronic density of fragment a would be an output 
variable of the dynamics of nucleus a. In general, how- 
ever, we find option (ii) to be more convenient: Fragment 
the static potential only, 



v(r) 



v a (r) 



(3) 



and define the time-dependent fragment potential v a (rt) 
by adding the total time-dependent potential u^(rt) to 
each of the v a (rt)'s; 

v a {rt) = v a (r) + v E (rt) . (4) 

Now define the many-electron fragment-a Hamiltonian 

as 

H a (t) = f + V ee + J d 3 r [v a (rt) + v p (rt)] h(r) . (5) 

The evolution of the state of this particular fragment is 
governed by the Liouville equation 



H a (t),t a (t) 



(6) 



The time-dependent electronic density of fragment a is 
given by n a (yt) = Tt{T a (t)n(r)} and the time- dependent 
partition potential v p (rt) of Eq. ^ is defined by requiring 
that the sum of fragment densities reproduce the total 
molecular density at all times: 



Nt 



n a (rt) = n(rt) 



(7) 



Just like traditional TDDFT is based on a one-to-one 
mapping between the Kohn-Sham potential v s (rt) and 
the electronic density n(rt), we now prove an analogous 
one-to-one mapping between n(rt) and v p (rt). The latter 
is therefore sharply defined by Eqs. ([T])-([7|. 

Theorem 1. For a given set of initial states {T a (to)}, 
the map between the density and the partition potential 
is invertible up to a time- dependent constant in the po- 
tential. 

Proof. The proof uses the Runge-Gross theorem pQ, and 
is analogous to it. Suppose there is a minimum integer 
k > such that 



d_ 

dt 



ml v' p (rt)-v p (rt)} 



= Constant 



< k 



t=t a I Constant m = k 



(8) 



Also assume that v p and v' p have the associated densi- 
ties {n Q } and {n' a } correspondingly. Suppose H a (t) and 
H' a {t) are the Hamiltonians of fragment a that corre- 
spond to v p and v' p , respectively. The key for the proof 
is the continuity equation 



dn a (rt) 
dt 



-V-ja(rt) 



(9) 



and the Liouville equation for the fragment current den- 
sities 



i^M=Tr 



dt 



r{t a (t) [j(r),ff a (t)]} . (10) 



Define 



Qk 



K(rt) - v p (rt)] 



(11) 



In virtue of the Runge-Gross theorem pQ and its gener- 
alization to ensembles [21], it is easy to show that 



d 



fc+2 



K(ri) - n a {rt)\ \ t=to 



Q t k+2 

Summing over all fragments gives 

Qk+2 



V-[n o (rt )V<(r)] . 

(12) 



Q t k+2 



[n'(rt)-n(rt)]| t=t0 =-V-[n(ri )V<(r)] . (13) 



Now we show that the right-hand side of this equation 
cannot be zero. Assume V • [7i(r£ )Vu^(r)] = and 
Vwj 7^ 0. Now invoke Green's identity to find 

J d 3 r^(r)V • (n(rto)Vti^(r)) = - J d 3 r n(rt )(Vw p k (r)) 2 



dS-n(rt )V(w p k ) 2 (r) = 



(14) 



If the total electronic density falls off enough to make 
the surface term negligible then Vw'(r) = 0, which is a 
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contradiction. Therefore the right-hand side of Eq. ( 13 ) the KS equations: 



cannot be zero. This leads to the conclusion that if v p 
and v p differ by more than a time-dependent constant 
then they cannot yield the same density in time. □ 



The above theorem shows that if {T a (to)} and v p (rt) 
are given, then one obtains a unique set of fragment den- 
sities {n a (rt)} and total density n(rt). The fragment 
density n a (rt) can be assumed to be non- interacting v- 
representable in time. Then we can associate a time- 
dependent Kohn-Sham potential w s , Q (ri) and initial state 
r s ,a(io) to describe the evolution of n a (rt) by means of 



idtfi a {rt) 



V s .a( rt ) 



where 



i(rt) = ^ n a( rt ) = ^2 f ia l^i«( r *)l' 



(15) 



(16) 



The fragments are allowed to have non-integer average 
numbers of electrons that are set by the initial state [17] . 
Since the Hamiltonian is particle-conserving, the occupa- 
tion numbers fi a remain fixed during the propagation. 

In analogy with PDFT, we define the xc potential by 
means of 



r«(to),r a , a (t )](ri) (*o)](ri) - v H [n a ](rt) - v p [n,{T a (t )}](rt) . (17) 



By comparing the fragment continuity equations for the 
interacting and non-interacting (Kohn-Sham) systems, 
we find that the above definition of the xc potential is 
consistent with (for example, see [22] ) 



V • < n a (rt)V 



3 n a (v't) 



d J r 



,(rt) 



V-[Q a (rt)-Q,, a (rt)] , 



(18) 



where the right-hand sides are hydrodynamical terms 
given by Q s ,a(rt) = -iTr{r s , Q (t)[j(r), f ]} and Q Q (rt) = 
-iTr{r a (t)[j(r),T+ V ee ]}- This indicates that the con- 
ventional xc potential of TDDFT and family of approxi- 
mations can be used for the fragments' TDKS equations, 
a direct consequence of van Leeuwen's theorem [23] , 

Furthermore, from the continuity equations for the to- 
tal current and fragment current densities, and from Eqs. 



(18), we find a formally exact relationship between the 



time-dependent partition potential and the total density: 



V ■ (n(rt)Vv p (rt)) 



d 2 n(rt) 
dt 2 



{v • Q sa [v P ](rt) 



-V ■ (n a [v p ](rt)Vv ata [v p ](rt))\ 

(19) 



where u s ,a[^p] = u hxc,<»K] + v a- In principle, evaluation 
of Eq. ( 19 ) at t — to yields a Sturm-Liouville linear differ- 
ential equation where v p (r,t — to) is the unknown vari- 
able. If we assume that the density is Taylor-expandable 
at t = to, then it is easy to show that consecutive differ- 
entiation of Eq. ( 19 ) and evaluation at t = to leads to 



a family of equations from which the Taylor coefficients 
of v p (rt) can be constructed in increasing order. This 
suggests that a given density is w p -representable as long 
as the conditions of the Sturm-Liouville theory are met. 

To illustrate our fragmentation approach, consider the 
simplest non-trivial model system consisting of a one- 
dimensional "electron" , two fragments, and an oscillating 
electric field of fixed frequency. For the static part of 
the external potential we choose a sum of soft-Coulomb 
potentials of equal strength Vq, a distance I apart: 



v{x) = V 



1 



+ 



y/(x + l/2)*+a y/(x-l/2)*+a 



(20) 

For the laser field we choose VE(x,t) — xE sm(ut), with 
E = 0.1, and u = 0.3. 

We partition the system by defining V\(x, t) = 
V /^/(x + l/2) 2 +a + v E (x,t) and v 2 (x,t) = 
V /^/(x-l/2) 2 + a + VE(x,t). The time-dependent 
fragment equations are, (for a = 1,2), 



i — <p a (xt) 



~2dx 2 



+ v a {xt) + v p {xt) <p a (xt) (21) 



The initial states of the fragments are obtained by solv- 
ing the ground-state PDFT equations as prescribed in 
Ref. [24]. This procedure generates the initial fragment 
Kohn-Sham orbitals needed to solve Eqs. [21] The dis- 
tance between the wells was chosen to allow for a signif- 
icant overlap between the initial fragments' densities. 
Even though the principle to construct v p (rt) is sim- 



ple, note that Eq. ( 19 1 can also be written as v p = J-v P 



where the operator T computes the right-hand-side of the 
equation, solves the differential equation, and finally out- 
puts v p (rt). One could employ this formula recursively, 
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Probability Densities 



Partition Potential v p {x, t) 




FIG. 1. Left: Fragment densities at times t = and t ~ T/4 
along with the total molecular density. Right: Exact partition 
potential v p (x,t). Parameters: I = 1, a = 1, Vo = — 1, 
E = 0.1, and cu = 0.3. A box of length 20 with 401 points 
was used. 



i.e 



vt: +1 = Fi 



We observed in our example that the 
term Q SQ becomes noisy even after short times if the sim- 
ulation box is discretized with large spatial steps. This 
noise is received by the partition potential during the 
propagation, and then it is received again by Q s a . This 
feedback process turns the algorithm unstable. The prob- 
lem is reminiscent to what occurs in traditional TDDFT 
when one wants to find the exchange-correlation poten- 
tial corresponding to a given density, even for only two 
electrons [25]. To solve this problem in TDDFT, Ref. [55] 
recently suggested an algorithm to control the feedback. 
They obtained encouraging results for a periodic system, 
but the methodology has not been tested for non-periodic 
systems. 

Instead, we found the exact time-dependent partition 
potential by using the following optimization procedure: 
The density and current density of the total system are 
found at each time step using the Crank-Nicolson prop- 
agator. (Other propagation methods may also be used.) 
A guess is made for the partition potential at the next 
unknown time and the fragment wave functions are prop- 
agated forward in time using this guess. (For small 
time steps the value of the partition potential at the 
previous time step works well.) The fragment densities 
({n a } , {ja}) are found using these fragment wave func- 
tions and added together to form an approximation to 
the total densities (n, j). The errors n CII = n — n cxact 
and jorr = j — Jcxact are computed and the residual 
norm (n orr /n CX act) + norm (j orr /j cxact ) is used in the L- 
BFGS-B optimizer [27], with the L 2 norm. The division 
by n oxact and jexact weights the error in the asymptotic 
regions to help increase the convergence rate, similar to 
the weighting used in [2"8] . 

The right panel of Fig. ([I]) displays the resulting par- 
tition potential. The left panel shows the total density, 
along the corresponding fragment densities at the initial 
time and at 1/ 4 th of a period. The importance of memory 
effects [29] is evident from Fig. ([2]), where the dashed- 
dotted lines labeled "Adiabatic" show the fragment den- 
sities obtained by solving the ground-state PDFT equa- 
tions for the instantaneous v(rt) at l/4 th of a period. 



Clearly, the correct partition potential is needed. Only 
when the electric field strength is reduced by a factor 
of 10 3 (keeping all other parameters fixed) does the adi- 
abatic partition potential produce a molecular density 
that is visibly indistinguishable from the exact molecular 
density at time t sa T/4. Interestingly, the approxima- 
tion v p (rt) ~ v p (rto) (labeled "Static" in Fig. ^) works 
qualitatively well for short times, certainly much better 
than the adiabatic approximation. The inset on the right 
panel of Fig. |2| shows how the static-w p approximation 
reproduces the correct dipole for short times. 
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FIG. 2. Left: Fragment densities at times t ~ T/4 along 
with the total molecular density (thick solid) calculated: with 
the exact Vp(xt) (solid), static v p (xt) = v p (xto) (dashed), and 
adiabatic v p (xt) = v p DFT [n(xt)] (dash-dotted). Right: Corre- 
sponding partition potentials at t ~ T/4. Parameters: 1 = 1, 
a = 1, Vo = —1, E — 0.1, and uj — 0.3. The inset compares 
the dipole moment obtained from the exact (solid) and static 
approximation to v p (xt) (dashed). The dotted line indicates 
the total (instantaneous) external potential. 



In practice, successful application of our approach to 
large systems will ultimately rely on the quality of ap- 
proximations to the time-dependent partition potential. 
The static approximation might be useful for short times. 
Furthermore, for problems whose physics is best de- 
scribed by invoking fragments (such as charge-transfer 
excitations), we believe that physically-meaningful ap- 
proximations for v p (rt) will be simpler to construct 
than approximations for the highly non-local exchange- 
correlation potential and kernel of TDDFT. Work along 
these lines, as well as on the linear-response formalism, 
is ongoing. 
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